Anomalous diffusion in supercooled liquids: a long-range localisation in particle 

trajectories 
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A statistical analysis oî the geometries oî particle trajectories in the supercooled liquid state is 
reported. We examine two structurally différent fragile glass-forming liquids simulated by molecular 
dynamics. In both liquids, the trajectories are found to exhibit a long-range localisation distinct 
from the short-range locahsation within the cage of nearest neighbours. This novel diffusion anomaly 
is interpreted as a resuit of the potential-energy landscape topography of fragile glass-formers where 
the local energy minima coalesce into metabasins - compact domains with low escape probability. 
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It is generally understood that the non-Arrhenius be- 
haviour and other dynamical anomalies observed in frag- 
ile ^ glass-forming liquids are caused by a particular to- 
pography of the potential energy-landscape (PEL) [3, Q ■ 
Having been cooled below the cross-over température Ta, 
the liquid stays within the basins of attraction of local 
PEL minima long enough for the velocities to decorre- 
late, and its structural relaxation is controUed by the 
inter-basin connectivity network. A distinctive feature 
of the PEL of a fragile liquid is that the basins are ge- 
ometrically organised into a two-scale pattern whereby 
compact sets of contiguous basins (metabasins) [3, [E] are 
separated by barriers with higher activation energy than 
those for elementary interbasin transitions. 

The impact of this PEL pattern upon the super- 
cooled liquid dynamics can be understood as foUows. 
In a multidimcnsional space, a boundary is close to al- 
most every point of the domain it confines Q. There- 
fore, almost every PEL minimum in a metabasin is 
separated by the metabasin's boundary from some of 
its Euclidean neighbour minima 0, [1]. The resulting 
réduction in the number of accessible relaxational de- 
grees of freedom slows down the relaxation and ren- 
ders elementary interbasin transitions increasingly col- 
lective yjlOj, in accordance with the Adam-Gibbs con- 
jecture [ll|. At a larger time scale, the system's confine- 
ment to a metabasin is expected to produce a localisa- 
tion effect in its configuration-space trajectory. It is as- 
sumed that the intra-metabasin dynamics correspond to 
/3-relaxation, whereas irréversible a-relaxation, and the 
respective decorrelation of the system's configuration- 
space trajectory, occur by virtue of the inter-metabasin 
transitions Q]- In contrast to collective intra-metabasin 
dynamics, relaxation beyond the metabasin boundaries 
is concluded to be mediated by uncorrelated movements 
of particles 12 1. 

The configuration-space trajectory of a liquid can be 
investigated by means of statistical analysis of its 3D- 



space projection - a trajectory of a single particle. Based 
on the above arguments, we can discriminate three dis- 
tinct régimes in the particle trajectory of a supercooled 
liquid (foUowing the initial ballistic régime), (i) Short- 
range localisation of a particle within the cage of its 
neighbours; (ii) Persistent diffusion due to the collective 
intra-metabasin dynamics 13]; (iii) Long-range localisa- 
tion produced by the the system's long-time confinement 
to a metabasin. Thèse anomalies can be detected as the 
respective déviations of Hausdorff measure from its value 
for the Brownian diffusion [3| . Régimes (i) and (ii) have 
been observed in molecular-dynamics simulations [Ï3| : 
the third one has not been reported so far. 

In this Letter, we analyse the statistical geometry of 
particle trajectories in the supercooled liquid state. We 
examine two structurally distinct fragile glass-forming 
liquids simulated by molecular dynamics. Evidence is 
presented for the long range localisation in the particles 
diffusion conjectured above as régime (iii). The observed 
effect can be regarded as direct évidence for the présence 
of metabasins in the PEL topography of fragile liquids. 
It also makes it possible to assess the spatial extent of 
a metabasin. We discuss the impact of the metabasin 
confinement on the relation between relaxation and dif- 
fusion. 

The MD simulations we report here explore two simple 
fragile glass-formers. One is the Z2 model [l5| demon- 
strating a pronounced tendency for icosahedral cluster- 
ing. This liquid was simulated at the number density 
p — 0.85. At that density, its estimated mode-coupling 
theory [3] critical température Te = 0.65, fragility in- 
dex B = 4.5 and Ta = 1.1- The liquid was simulated at 
two températures: T = 1.2 and T = 0.7, above and be- 
low Ta, respectively. The other glass-forming liquid ex- 
plored here was the binary Lennard-Jones System (BLJ) 



commonly used in glass studies [17 . 
action potential was truncated at r 
Stroddard-Ford quadratic cutoff [ic 



IJ, |19|. The inter- 
= 3.0036 using the 
|. The BLJ hquid 
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was simulated in a supercooled state at p = 1.19 and 
T = 0.451. In this model, we analysed the trajectories of 
thc smaller (more mobile) particles. Both models were 
simulated in the NVE ensemble with N = 16000. 

Traditionally, statistics of diffusing particles is pre- 
sented in terms of the distribution of partiale displace- 
ments at time t, Gs{r,t) [lïj. In the case of Brow- 
nian diffusion, this is a Gaussian, and its second mo- 
ment, the mean-square displacement, dépends on t as 
(r(t)^} = 6Dt, D being the diffusion coefficient. 

Alternatively, the particles' diffusion can be described 
in terms of the first-passage time distribution P{e, t) [2^ . 
P{e,t)dt is the probability that a particle's first crossing 
of the boundary of a sphère of radius e centered at the 
particle's initial position happens within the time interval 
between t and t + dt. For Brownian diffusion, the first 
moment of this distribution, the mean first-passage time 
Te is equal to = e'^/{6D) 22]. 

In this study, we analyse the statistical geometry of 
particle trajectories. It is convenient for this purpose to 
use the distribution of first passage trajectory lengths 
(FPL) P{e,L), rather than first passage times. Thèse 
two distributions are équivalent for the description of liq- 
uid diffusion (except for very short times) since, in dense 
liquids, velocities decorrelate much faster than trajecto- 
ries. The length of a particle trajectory is defined as the 
time intégral of its instantaneous speed. With good accu- 
racy, the trajectory length within time t is L = vt, where 
V = ^/8T/^T is the mean particle speed [i^. We dénote 
the first moment of the FPL as L^. For the trajectory of 
a Brownian particle we have 
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In this study we explore corrélations in the system's 
configuration-space trajectory and use the decay of the 
corrélations as an indicator of the system's approach to 
ergodic equilibrium. The présence of corrélation effects 
in a particle trajectory can be detected from its scaling 
behaviour. The latter can be quantitatively analyzed in 
terms of the trajectory's fractal dimensionality (Haus- 
dorff measure) [1J|. The Hausdorff measure a is defined 
as the logarithmic derivative of : 
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As it foUows from équations ([T]) and ([2]), a(e) — 2 
for a particle trajectory in thc largc-e limit of Brownian 
diffusion. 

Within the intermediate range of e, the déviation of 
a{e) from the indicated Brownian value can be used 
to discriminate the spécifie diffusion régimes discussed 
above. If a{e) > 2, consécutive particle displacements, 
defined by the scale of e, correlate negatively (the dif- 
fusion is localised). If a{e) < 2, consécutive particle 



displacements correlate positively; this is a signature of 
persistent diffusion. 

Fig. 1 shows a{e) for thc three liquid states simu- 
lated in this study. A common feature in ail cases is a 
pronounced peak located at around e = 0.3 where a{e) 
exceeds the Brownian value a{e) = 2. The peak, greatly 
enhanced by cooling below T^, indicates négative cor- 
relations in the particle trajectory within the respective 
range of e. It apparently represents a localisation effect 
induced by the particle's confinement to the cage of near- 
est neighbours. This can be seen as a 3D projection of 
the configuration-space trajectory localisation within a 
PEL minimum basin. 

For the Z2 liquid at T = 1.2, representing normal liquid 
dynamics, the described short-range localisation régime 
is immediately foUowed at larger e by convergence to the 
asymptotic Brownian behaviour. In the case of the super- 
cooled liquid dynamics, both for the Z2 liquid at T = 0.7 
and the BLJ liquid, a broad minimum is observed within 
0.8 < e < 3 where a{e) < 2. This effect, indicating 
positive corrélations in a particle trajectory, is character- 
istic of persistent diffusion (régime (ii) discussed above). 
A similar behaviour has already been reported for the 
BLJ liquid, and concluded to be related to the collective 
particle motions [ï^. 

At still larger values of e we observe, both for the Z2 
liquid at T = 0.7 and the BLJ liquid, another anomaly 
in the behaviour of a{e). Before eventually converging 
to the Brownian limit, a{e) in thèse supercooled liquids 
exhibits a second maximum located within 3 < e < 10 
where a{s) > 2. This feature, shown in the inset in 
Fig. 1, indicates the présence of a long-range localisation 
in the particle diffusion that we conjectured above as 
régime (iii). This interesting novel effect is the central 
observation of the présent study, and we shall discuss it 
in more détail. 

Besides Brownian diffusion, another essential criterion 
of the ergodic behaviour of a liquid is homogeneity of 
the diffusion process. In terms of the statistics of the 
particles' trajectories, this can be measured by compar- 
ing the variance of the FPL distribution: Var[P{L, e)] = 
(L(e)^)/I/g — 1 with its asymptotic value of 2/5 expected 
for Brownian diffusion [22{ . The évolution of this quan- 
tity with respect to e for ail three simulated liquid states 
is shown in Fig. 2. Comparing thèse data with those 
in Fig. 1 makes it possible to conclude that the anoma- 
lies in the fractal dimensionality observed in the super- 
cooled liquid state, both for Z2 and BLJ, are accompa- 
nied by anomalously high levels of variance of the FPL 
distribution. Excessive variance is observed in both the 
localisation and the persistent diffusion domains. This 
observation is consistent with the point of view that the 
anomalous diffusion in the supercooled liquid state is as- 
sociated with dynamical heterogeneity. 

As we argued above, the localisation régime that we 
observed as a large-e maximum in a(e) for the super- 
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FIG. 1: Apparent fractal dimension of a particle trajectory. 
Solid black line and dashed red line: Z2 liquid for T = 0.7 and 
T = 1.2, respectively. Dash-dotted green line, the BLJ liquid. 
The errer bar indicate the estimated standard déviation for 
the BLJ data. 
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FIG. 2: The variance of the FPL distribution. Sohd black 
line and dashed red line: Z2 liquid for T = 0.7 and T = L2, 
respectively. Dash-dotted green line, the BLJ liquid. 



cooled liquid state in Fig. 1 can be viewed as a resuit 
of the system's confinement to a configuration-space do- 
main that we associate with a metabasin. In this way, 
the effect in question is a direct observation of the con- 
straints imposed by the metabasin structure of the PEL 
upon the liquid diffusion, and its location in tcrms of e 
makes it possible to conclude the extent of that structure. 

It is worth noting that the observed long-range maxi- 
mum of a(e) makes the diffusion at shorter scales faster 
as compared with the asymptotic Brownian diffusion. To 
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FIG. 3: Ratio of the e-dependent diffusion rate (Eq. 3) 
to the asymptotic Brownian diffusion rate. Solid black line 
and dashed red line: Z2 liquid for T = 0.7 and T = 1.2, 
respectively. Dash-dotted green line, the BLJ liquid. 
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FIG. 4: FPL distribution for the particle trajectories in Z2 
liquid for e = 2.424. L is scaled by ve^ / (6D), the e-dependent 
mean first passage trajectory length for Brownian diffusion 
(Eq. 1), D being the asymptotic diffusion rate for the respec- 
tive T. Solid line, T = 0.7; dashed red line, T = 1.2. The 
inset depicts, in log-linear scale, the large-L behaviour of the 
FPL. 

demonstrate this effect, we introduce the e-dependent dif- 
fusion rate using Eq. 1: 

D, = vey{6L,) (3) 

Fig. 3 displays the ratio of D^/D, where D is the 
Brownian diffusion rate derived from the asymptotic be- 
haviour of the mean-square displacement in the limit of 
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large time (in that limit, and D converge). The gên- 
erai behaviour oî D^/D is related to that of a{e), Fig. 1, 
by Eq. 2. The long-range maxima of a{e) observed for 
both supercooled liquids, although appear subtle, obvi- 
ously resuit in a measurable réduction of D^/ D within 
the same range of e, as can be seen in Fig. 3. 

To obtain further insight into this effect, we compare in 
Fig. 4 two FPL distributions for the Z2 liquid, at T = 0.7 
and at T = 1.2, in both cases for e — 2.424. As can be 
seen in Fig. 3, at T = 0.7 that value of e represents 
the diffusion enhancement régime, whereas for T = 1.2 
the diffusion appears to attain Brownian behaviour. In 
order to compare the distributions on equal ground, we 
scaled L by ve'^/{6D), the e-dependent mean first passage 
trajectory length for the Brownian diffusion (Eq. 1), D 
being the asymptotic diffusion rate for the respective T. 
The large-L tails of the distributions (shown in the inset) 
are practically indistinguishable, and their asymptotic 
behaviour is apparcntly exponential. The différence in 
the first moments of thèse two distributions arises from a 
significant extra contribution from shorter-L trajectories 
in the case of T = 0.7 indicating an excessive présence 
in the supercooled liquid of highly mobile particles. This 
observation is apparently consistent with both the diffu- 
sion enhancement as observed in Fig. 3, and the excess 
variance in the FPL distribution within the respective 
range of e, Fig. 2. 

A conceptually interesting aspect of thèse results con- 
cerns the relation between the diffusion and relaxation 
[il ]. How far does a particle of a liquid diffuse before 
the ergodic equilibrium has been attained? A necessary 
condition for that is the onset of an equilibrium (Brown- 
ian) diffusion régime in the particle trajectories. We see 
in Figs. 1 and 3 that, in a normal liquid state above 
Ta, Brownian diffusion is reached as soon as a parti- 
cle leaves the cage of its neighbours. The appearance 
of two longer-range anomalous diffusion régimes in the 
supercooled liquid state, which we attribute to the con- 
finement to a metabasin, extends the spatial scale of 
the non-equilibrium diffusion process, thereby delaying 
the relaxation relative to the diffusion. In this way, the 
latter anomaly can also be regarded as an effect of the 
metabasin topography of the PEL of fragile liquids. 

In summary, we find that two physically and struc- 
turally différent supercooled fragile glass-forming liquids 
demonstrate a long-range localisation in the particle tra- 
jectories. We link this novel effect to the existence 
of mctabasins in the fragile liquids' PEL topographies. 
Based on the results of this study, and our arguments pre- 
sented above, the gênerai impact of the PEL mctabasins 
upon the supercooled liquid dynamics can be described 
as foUows. First, constraints imposcd by the metabasin 
boundary reduce the number of accessible degrees of 
freedom for the intra-metabasin transitions between the 
basins of the PEL minima. This induces, within the re- 
spective range of distances, collective dynamics and, con- 



sequently, persistent diffusion with a(e) < 2. Second, on 
a larger spatial scale, the system's (time-limited) con- 
finement to a compact configuration-space domain de- 
fined as a metabasin results in a long-range localisation 
of the particle trajectories. In this régime, a{e) > 2. 
We also show that thèse anomalies effectively give an 
intermediate-range diffusion rate that is faster than the 
asymptotic Brownian diffusion. Finally, the existence of 
the long-range anomalies in a particle trajectory, delay- 
ing its approach to the equilibrium Brownian diffusion, 
also delays the liquid's relaxation relative to diffusion. 
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